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TITIJE OF THE INVENTION 
METHOD OF FAULT DETECTION FOR 
MATERIAL PROCESS SYSTEM 

CROSS-REFERENCE TO RELATED APPLICATIONS 
[0001] This application claims priority and is related to United States application serial no. 
60/343,175, filed on December 31, 2001, the contents of which are herein incorporated by 
reference. This application is related to co-pending PCT application serial no. 
PCT/US02/XXXXX, filed on even date herewith, Attorney Docket No. 216951WO, the 
contents of which are herein incorporated by reference. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0002] The present invention relates to material processing and more particularly to a method 
for fault detection for material process system. 

Description of Related Art 

[0003] One area of material processing in the semiconductor industry which presents 
formidable challenges is, for example, the manufacture of integrated circuits (ICs). Demands 
for increasing the speed of ICs in general, and memory devices in particular, force 
semiconductor manufacturers to make devices smaller and smaller on the wafer surface. And 
conversely, while shrinking device sizes on the substrate is incurred, the number of devices 
fabricated on a single substrate is dramatically increased with further expansion of the 
substrate diameter (or processing real estate) from 200 mm to 300 mm and greater. Both the 
reduction in feature size, which places greater emphasis on critical dimensions (CD), and the 
increase of substrate size lead to even greater requirements on material processing uniformity 
to maximize the yield of superior devices. 
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[0004] Typically, during materials processing, one method to facilitate the addition and 
removal of material films when fabricating composite material structures includes, for 
example, the use of plasma. For example, in semiconductor processing, a (dry) plasma etch 
process is utilized to remove or etch material along fine lines or within vias or contacts 
patterned on a silicon substrate. 

[0005] During, for example, material processing in IC fabrication, shrinking critical feature 
sizes, increasing substrate sizes and escalating numbers and complexities of processes lead to 
the necessity to control the material processing uniformity throughout the lifetime of a 
process and from process-to-process. The lack of uniformity in simply one process 
measurable generally requires the sacrifice of other important process parameters, at least, 
somewhere during the process. In material processing, the lack of process uniformity can, for 
example, cause a costly reduction in the yield of superior devices. 

[0006] Attempts to design material processing hardware either to produce uniform processing 
properties or correct for known non-uniformities are further complicated by the expansive set 
of independent parameters, the complexity of these material processing devices, and simply 
the exorbitant cost and lack of robustness of such material processing devices. Furthermore, 
for conventional material processing devices, the number of externally, controllable 
parameters are severely limited to only a few known, adjustable parameters. Therefore, it is 
essential that the inter-relations between all externally controllable parameters and 
measurable process parameters are derived and made useful throughout the lifetime of a 
process and from process-to-process. 
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SUMMARY CW TTTE INVENTION 
[0007] The present invention provides for a method of material processing, wherein the 

method comprises utilizing a material processing system, the material processing system 

comprising a process chamber, a device for measuring and adjusting at least one controllable 

process parameter, and a device for measuring at least one process performance parameter. 

[0008] The present invention provides a method of characterizing a process, the 

characterizing comprising identifying a signature of the process wherein the signature 

comprises at least one spatial component; optimizing the process, the optimizing comprising 

identifying a reference signature; comparing the signature of the process with the reference 

signature for the process, wherein the comparing comprises determining a difference 

signature, and determining a process fault by comparing the difference signature with a 

threshold, wherein the process fault occurs when the threshold is exceeded. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0009] These and other objects and advantages of the invention will become more apparent 

and more readily appreciated from the following detailed description of the exemplary 
embodiments of the invention taken in conjunction with the accompanying drawings, where: 
[0010] FIG. 1 shows a material processing system according to a preferred embodiment of 
the present invention; 

[0011] FIG. 2 shows a material processing system according to an alternate embodiment of 
the present invention; 

[0012] FIG. 3 shows a material processing system according to another embodiment of the 
present invention; 

[0013] FIG. 4 shows a material processing system according to another embodiment of the 
present invention; 
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[0014] FIG. 5 shows a material processing system according to an additional embodiment of 
the present invention; 

[0015] FIG. 6A presents a data scan of a first etch rate profile; 

[0016] FIG. 6B presents a spectrum of spatial components for the data scan of FIG. 6A; 
[0017] FIG. 7A presents a data scan of a second etch rate profile; 
[0018] FIG. 7B presents a spectrum of spatial components for the data scan of FIG. 7A; 
[0019] FIG. 8A presents a comparison of the spectrum of spatial components resulting from 
an increase in process pressure; 

[0020] FIG. 8B presents the difference spectrum for the data of FIG. 8A; 

[0021] FIG. 9A presents a comparison of the spectrum of spatial components resulting from a 

decrease in RF power; 

[0022] FIG. 9B presents the difference spectrum for the data of FIG. 9A; 

[0023] FIG. 9C presents the difference spectrum for an increase in RF power; 

[0024] FIG. 10A shows an exemplary spectrum of spatial components for a non-uniform etch 

rate; 

[0025] FIG. 10B shows an exemplary spectrum of spatial components for a uniform etch 
rate; 

[0026] FIG 11 presents an exemplary table of variations in spatial components provided 
changes in controllable process parameters; 

[0027] FIG. 12 presents an exemplary plot of the cumulative sum of squares and cumulative 
sum of variations to the sum of squares for three principal components; 
[0028] FIG. 13A presents the scores corresponding to each spatial component in t(l), t(2) 
space provided the exemplary data of FIG. 11; 

[0029] FIG. 13B presents the loadings for each variable in p(l), p(2) space provided the 
exemplary data of FIG. 11; 
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[0030] FIG. 14A presents the scores corresponding to each spatial component in t(l), t(3) 
space provided the exemplary data of FIG. 11; 

[0031] FIG. 14B presents the loadings for each variable in p(l), p(3) space provided the 
exemplary data of FIG. 11; 

[0032] FIG. 15 presents an exemplary table summarizing data presented in FIGs. 13A,B and 
14A,B; 

[0033] FIG. 16A presents a table of spatial components for a reduced set of the data of the 
table presented in FIG. 1 1 ; 

[0034] FIG. 16B presents a spectrum of spatial components according to the data of FIGs. 
6A,B, and a spectrum of spatial components according to the data of FIG. 16A; 
[0035] FIG. 16C presents a difference spectrum obtained from the spectra of FIG. 16B; 
[0036] FIG. 17 shows a data scan of a first etch profile according to the data of FIGs. 6A,B, 
and a data scan of a second etch profile according to the data of FIG. 16C; 
[0037] FIG. 18 A presents a flow diagram of a method according to the present invention; 
[0038] FIG. 18B presents a flow diagram of an additional method according to the present 
invention; 

[0039] FIG. 18C presents a flow diagram of an additional method according to the present 
invention; and 

[0040] FIG. 18D presents a flow diagram of an additional method according to the present 
invention. 

DRTATT FD DESCRIPTION OF EXEMPLARY EMBODIMENTS 
[0041] According to an embodiment of the present invention, a material processing system 1 

is depicted in FIG. 1 comprising a process chamber 10, a device for measuring and adjusting 

at least one controllable process parameter 12, a device for measuring at least process 

performance parameter 14, and a controller 55. The controller 55 is coupled to the device for 
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measuring and adjusting at least one controllable process parameter 12 and the device for 
measuring at least one process performance parameter. Moreover, the controller 55 is 
capable of executing the method of performing a process to be described. 
[0042] In the illustrated embodiment, material processing system 1, depicted in FIG. 1, 
utilizes a plasma for material processing. Desirably, material processing system 1 comprises 
an etch chamber. Alternately, material processing system 1 comprises a photoresist coating 
chamber such as, for example, a photoresist spin coating system. In another embodiment, 
material processing system 1 comprises a photoresist patterning chamber such as, for 
example, a ultraviolet (UV) lithography system. In another embodiment, material processing 
system 1 comprises a dielectric coating chamber such as, for example, a spin-on-glass (SOG) 
or spin-on-dielectric (SOD) system. In another embodiment, material processing system 1 
comprises a deposition chamber such as, for example, a chemical vapor deposition (CVD) 
system or a physical vapor deposition (PVD) system. In an additional embodiment, material 
processing system 1 comprises a rapid thermal processing (RTP) chamber such as, for 
example, a RTP system for thermal annealing. In another embodiment, material processing 
system 1 comprises a batch diffusion furnace. 

[0043] According to the illustrated embodiment of the present invention depicted in FIG. 2, 
material processing system 1 can comprise process chamber 10, substrate holder 20, upon 
which a substrate 25 to be processed is affixed, gas injection system 40, and vacuum 
pumping system 50. Substrate 25 can be, for example, a semiconductor substrate, a wafer or 
a liquid crystal display. Process chamber 10 can be, for example, configured to facilitate the 
generation of plasma in processing region 45 adjacent a surface of substrate 25, wherein 
plasma is formed via collisions between heated electrons and an ionizable gas. An ionizable 
gas or mixture of gases is introduced via gas injection system 40 and the process pressure is 
adjusted. For example, a control mechanism (not shown) can be used to throttle the vacuum 
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pumping system 50. Desirably, plasma is utilized to create materials specific to a pre- 
determined materials process, and to aid either the deposition of material to substrate 25 or 
the removal of material from the exposed surfaces of substrate 25, 

[0044] Substrate 25 can be, for example, transferred into and out of chamber 10 through a 
slot valve (not shown) and chamber feed-through (not shown) via robotic substrate transfer 
system where it is received by substrate lift pins (not shown) housed within substrate holder 
20 and mechanically translated by devices housed therein. Once substrate 25 is received 
from substrate transfer system, it is lowered to an upper surface of substrate holder 20. 
[0045] Desirably, the substrate 25 can be, for example, affixed to the substrate holder 20 via 
an electrostatic clamping system 28. Furthermore, substrate holder 20 can further include a 
cooling system including a re-circulating coolant flow that receives heat from substrate 
holder 20 and transfers heat to a heat exchanger system (not shown), or when heating, 
transfers heat from the heat exchanger system. The heating/cooling system further comprises 
a device 27 for monitoring the substrate 25 and/or substrate holder 20 temperature. The 
device 27 can be, for example, a thermocouple (e.g. K-type thermocouple), pyrometer, or 
optical thermometer. Moreover, gas can be delivered to the back-side of the substrate via a 
backside gas system 26 to improve the gas-gap thermal conductance between substrate 25 
and substrate holder 20. Such a system can be utilized when temperature control of the 
substrate is required at elevated or reduced temperatures. For example, temperature control 
of the substrate can be useful at temperatures in excess of the steady-state temperature 
achieved due to a balance of the heat flux delivered to the substrate 25 from the plasma and 
the heat flux removed from substrate 25 by conduction to the substrate holder 20. In other 
embodiments, heating elements, such as resistive heating elements, or thermo-electric 
heaters/coolers can be included. 
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[0046] In the illustrated embodiment, shown in FIG. 2, substrate holder 20 can, for example, 
further serve as an electrode through which RF power is coupled to plasma in processing 
region 45. For example, substrate holder 20 is electrically biased at a RF voltage via the 
transmission of RF power from RF generator 30 through impedance match network 32 to 
substrate holder 20, The RF bias can serve to heat electrons and, thereby, form and maintain 
plasma. In this configuration, the system can operate as a reactive ion etch (RE) reactor, 
wherein the chamber and upper gas injection electrode serve as ground surfaces. A typical 
frequency for the RF bias can range from 1 MHz to 100 MHz (e.g., 13.56 MHz). RF systems 
for plasma processing are well known to those skilled in the art. 

[0047] Alternately, RF power is applied to the substrate holder electrode at multiple 
frequencies. Furthermore, impedance match network 32 serves to maximize the transfer of 
RF power to plasma in processing chamber 10 by minimizing the reflected power. Match 
network topologies (e.g. L-type, Tt-type, T-type, etc.) and automatic control methods are well 
known to those skilled in the art. 

[0048] With continuing reference to FIG. 2, process gas 42 can be, for example, introduced 
to processing region 45 through gas injection system 40. Process gas 42 can, for example, 
comprise a mixture of gases such as argon, CF 4 and 0 2 , or argon, C4F8 and 0 2 for oxide etch 
applications. Gas injection system 40 can comprise a showerhead, wherein process gas 42 is 
supplied from a gas delivery system (not shown) to the processing region 45 through a gas 
injection plenum (not shown), a series of baffle plates (not shown) and a multi-orifice 
showerhead gas injection plate (not shown). Gas injection systems are well known to those 
of skill in the art. 

[0049] Vacuum pump system 50 can, for example, include a turbo-molecular vacuum pump 
(TMP) capable of a pumping speed up to 5000 liters per second (and greater) and a gate valve 
for throttling the chamber pressure. In conventional plasma processing devices utilized for 
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dry plasma etch, a 1000 to 3000 liter per second TMP is employed. TMPs are useful for low 
pressure processing, typically less than 50 mTorr. At higher pressures, the TMP pumping 
speed falls off dramatically. For high pressure processing (i.e. greater than 100 mTorr), a 
mechanical booster pump and dry roughing pump can be used. Furthermore, a device for 
monitoring chamber pressure 52 is coupled to the chamber 10. The pressure measuring 
device 52 can be, for example, a Type 628B Baratron absolute capacitance manometer 
commercially available from MKS Instruments, Inc. (Andover, MA). 

[0050] Material processing system 1 further comprises a metrology tool 100 to measure 
process performance parameters such as, for example for etch systems, an etch rate, an etch 
selectivity (i.e. ratio of etch rate of one material to etch rate of a second material), an etch 
uniformity, a feature profile angle, a critical dimension, etc. The metrology tool 100 can be 
either an in-situ or ex-situ device. For an in-situ device, the metrology tool 100 can be, for 
example, a scatterometer, incorporating beam profile ellipsometry and beam profile 
reflectometry, commercially available from Therma-Wave, Inc. (1250 Reliance Way, 
Fremont, CA 94539) which is positioned within the transfer chamber (not shown) to analyze 
substrates 25 transferred into and out of process chamber 10. For an ex-situ device, the 
metrology tool 100 can be, for example, a scanning electron microscope (SEM) wherein 
substrates have been cleaved and features are illuminated to determine the above performance 
parameters. The latter approach is well known to those skilled in the art of substrate 
inspection. The metrology tool is further coupled to controller 55 to provide controller 55 
with spatially resolved measurements of the process performance parameters. 
[0051] Controller 55 comprises a microprocessor, memory, and a digital I/O port capable of 
generating control voltages sufficient to communicate and activate inputs to material 
processing system 1 as well as monitor outputs from material processing system 1. 
Moreover, controller 55 is coupled to and exchanges information with RF generator 30, 
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impedance match network 32, gas injection system 40, vacuum pump system 50, pressure 
measuring device 52, backside gas delivery system 26, substrate/substrate holder temperature 
measurement system 27, electrostatic clamping system 28, and metrology tool 100. A 
program stored in the memory is utilized to activate the inputs to the aforementioned 
components of a material processing system 1 according to a stored process recipe. One 
example of controller 55 is a DELL PRECISION WORKSTATION 610™ , available from 
Dell Corporation, Dallas, Texas. 

[0052] In the illustrated embodiment, shown in FIG. 3, the material processing system 1 can, 
for example, further comprise either a mechanically or electrically rotating DC magnetic field 
system 60, in order to potentially increase plasma density and/or improve plasma processing 
uniformity, in addition to those components described with reference to FIGs. 1 and 2. 
Moreover, controller 55 is coupled to rotating magnetic field system 60 in order to regulate 
the speed of rotation and field strength. The design and implementation of a rotating 
magnetic field is well known to those skilled in the art. 

[0053] In the illustrated embodiment, shown in FIG. 4, the material processing system 1 of 
FIGs. 1 and 2 can, for example, further comprise an upper electrode 70 to which RF power 
can be coupled from RF generator 72 through impedance match network 74. A typical 
frequency for the application of RF power to the upper electrode can range from 10 MHz to 
200 MHz (e.g., 60 MHz). Additionally, a typical frequency for the application of power to 
the lower electrode can range from 0.1 MHz to 30 MHz (e.g., 2 MHz). Moreover, controller 
55 is coupled to RF generator 72 and impedance match network 74 in order to control the 
application of RF power to upper electrode 70. The design and implementation of an upper 
electrode is well known to those skilled in the art. 

[0054] In the illustrated embodiment, shown in FIG. 5, the material processing system of 
FIG. 1 can, for example, further comprise an inductive coil 80 to which RF power is coupled 
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via RF generator 82 through impedance match network 84. RF power is inductively coupled 
from inductive coil 80 through dielectric window (not shown) to plasma processing region 
45. A typical frequency for the application of RF power to the inductive coil 80 can range 
from 10 MHz to 100 MHz (e.g., 13.56 MHz). Similarly, a typical frequency for the 
application of power to the chuck electrode can range from 0.1 MHz to 30 MHz (e.g., 13.56 
MHz). In addition, a slotted Faraday shield (not shown) can be employed to reduce 
capacitive coupling between the inductive coil 80 and plasma. Alternately, coil 80 can be 
positioned above chamber 10 as a spiral-like coil such as in a transformer coupled plasma 
(TCP) source. Moreover, controller 55 is coupled to RF generator 82 and impedance match 
network 84 in order to control the application of power to inductive coil 80. The design and 
implementation of an inductively coupled plasma (ICP) source and a transformer coupled 
plasma (TCP) source are well known to those skilled in the art. 

[0055] Alternately, the plasma can be formed using electron cyclotron resonance (ECR). In 
yet another embodiment, the plasma is formed from the launching of a Helicon wave. In yet 
another embodiment, the plasma is formed from a propagating surface wave. Each plasma 
source described above is well known to those skilled in the art. 

[0056] Referring now to FIGs. 1 through 5, substrates 25 are processed in process chamber 
10 and some process performance parameters can be measured utilizing, for example, the 
metrology tool 100. Desirably, process performance parameters can include, for instance, 
etch rate, deposition rate, etch selectivity (ratio of the rate at which a first material is etched 
to the rate at which a second material is etched), an etch critical dimension (e.g. length or 
width of feature), an etch feature anisotropy (e.g. etch feature sidewall profile), a film 
property (e.g. film stress, porosity, etc.), a plasma density (obtained, for example, from a 
Langmuir probe), an ion energy (obtained, for example, from an ion energy spectrum 
analyzer), a concentration of a chemical specie (obtained, for example, from optical emission 
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spectroscopy), a temperature, a pressure, a mask (e.g. photoresist) film thickness, a mask (e.g. 
photoresist) pattern critical dimension, etc. For example, FIG. 6A presents a substrate scan of 
the etch rate (Angstroms/minute, A/min) as a function of position (millimeters, mm) on a first 
substrate 25, where a position of zero (0) corresponds to the center of substrate 25 and a 
position of plus or minus (±) 100 corresponds to diametrically opposite edges of a, for 
example, (200 mm) substrate 25. Similarly, FIG. 7A presents a substrate scan of the etch rate 
versus substrate position for a second substrate 25. 

[0057] In FIGs. 6A and 7A, thirty-two (32) samples are taken along a full radial scan (edge- 
to-edge) of the substrate diameter; however, in general, the number of samples can be 
arbitrary, e.g. N samples where N^2. The time T required to input the data at a sampling 
rate R can be expressed as T=N/R; i.e. T=N/R=(32 samples)/(1000 samples/second)=0.032 
seconds (for sampling 32 points across a substrate at 1 kHz). For a data scan of period T, the 
primary spatial component is f=l/T and the highest spatial component must satisfy the 
Nyquist critical frequency of fmax^l^A, where A=T/N. Therefore, in the above example, 
f=l/T=R/N=31.25 Hz and W=l/2A==R/2==500 Hz. 

[0058] In general, a scan of data, as described above, can be manipulated into spectral space 
and be represented by a set of orthogonal components. For example, if the samples are 
equally spaced in time (or space) and the scan is assumed to be periodic, then the data scan is 
directly amenable to the application of a discrete Fourier transform of the data to convert the 
data from physical space to Fourier (spectral) space. Moreover, if the samples are unequally 
spaced in time (or space), there exist methods of treating the data. These methods are known 
to those skilled in the art of data processing. When using a Fourier series representation of 
the data, the spatial components can be, for example, Fourier harmonics. Moreover, if the 
sampling period T is relatively small (small relative to the change of the data scan in time; 
applicable only for in-situ monitoring during substrate processing), then the Fourier spectrum 
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can be regarded as a wavenumber spectrum and the minimum and maximum spatial 
components can be referred to as the minimum and maximum wavenumbers (or maximum 
and minimum wavelengths, respectively). 

[0059] FIG. 6B presents the amplitude of each spatial component (i.e. f n =n/NA, n=l,N/2) for 
the data scan shown in FIG. 6A. And similarly, FIG. 7B presents the amplitude of each 
spatial component (i.e. f n =n/NA, n=l,N/2) for the data scan shown in FIG. 7 A. In general, 
the following observations can be made: (1) the primary spatial component (f x ) has the largest 
magnitude and represents contributions from each point in the data scan (therefore, all points 
are interdependent; longest wavelength); and (2) the highest spatial component (f m ) has 
typically the smallest magnitude and it represents each point in the data scan separately 
(therefore, all points are independent of one another, smallest wavelength). Additionally, 
subtle changes in the etch rate profile (i.e. FIG. 6A versus FIG. 7A) can have a significant 
effect on the signature described by the spatial components in spectral space (i.e. FIG. 6B 
versus 7B). 

[0060] Therefore, changes in the signature (spectrum) of spatial components can indicate 
whether the process variations leading to the observed spectral shifts are occurring globally 
over the substrate or locally over the substrate. In summary, changes in the amplitudes of the 
lower order spatial components (i.e. f u f 2 , f 3 , ..-) reflect global variations of processing 
parameters above substrate 25, and changes in the amplitudes of the higher order spatial 
components (i.e. ... ,f N /2- 2 , fww. fiw) reflect local variations of processing parameters above 
substrate 25. 

[0061] For example, a change in the pressure or RF power (e.g. an increase in the processing 
pressure or decrease in the RF power) is expected to have a global effect on the signature of 
spatial components and, hence, affect primarily the lower order components. FIG. 8A 
presents an example of raising the chamber pressure and its effect on the signature of spatial 
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components, where FIG. 8B presents the respective difference signature. Similarly, FIG. 9A 
presents an example of reducing the RF power and its effect on the signature of spatial 
components, where FIG. 9B presents the respective difference signature (for reducing the RF 
power) and FIG. 9C presents the corresponding difference signature for increasing the RF 
power. Each difference signature provides a different spatial characteristic (i.e. 'Yingerprint") 
for each type of process change (i.e. increase or decrease in process pressure, increase or 
decrease in RF power, increase or decrease in mass flow rate of process gas, etc.). 
[0062] Since each process performed in material processing system 1 can be characterized by 
a signature of its spatial components, then one can evaluate the effect of process uniformity 
on the signature of spatial components. FIG. 10A presents a spatial signature for a non- 
uniform process and FIG. 10B presents a spatial signature for a uniform process. Clearly, the 
uniformity of a process can be directly correlated with an overall reduction in the magnitudes 
of each spatial component. 

[0063] Since there exists a relationship between controllable process parameters and spectra 
of spatial components obtained from a scan of, for instance, the etch rate across the substrate, 
then it is conceivable that spatial component differences can be subjected to linear 
superposition, i.e added and subtracted, to minimize the magnitudes of all spatial components 
and, therefore, produce a uniform process. A method of establishing a correlation between 
changes in controllable process parameters and spatial components utilizing multivariate 
analysis is now described to determine the right combination of variables to produce a 
uniform process. 

[0064] The table provided in FIG. 1 1 presents the relative change in amplitude of each spatial 
component through the first sixteen (16) components for twelve (12) variations in 
controllable process parameters. The controllable process parameters include, for example: 
(1) increase in process pressure, (2) decrease in process pressure, (3) increase in (Helium) 
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backside gas pressure, (4) decrease in (Helium) backside gas pressure, (5) increase in CF 4 
partial pressure, (6) decrease in CF 4 partial pressure, (7) increase in RF power, (8) decrease in 
RF power, (9) increase in substrate temperature, (10) decrease in substrate temperature, (11) 
the use of a 12 mm focus ring, and (12) the use of a 20 mm focus ring (instead of the nominal 
16 mm focus ring). Each of the above exemplary controllable process parameters are 
measurable and adjustable with reference to FIGs. 1 through 5. The process pressure can be 
adjusted and monitored during process using either changes in, for example, the gate valve 
setting or the total process gas mass flow rate, in concert with a pressure measuring device 
52. The forward and reflected RF power can be adjusted and monitored using commands to 
the RF generator 30 (FIG. 2), the match network 32 (FIG. 2), a dual directional coupler (not 
shown) and power meters (not shown). The CF 4 partial pressures can be adjusted and 
monitored using a mass flow controller to regulate the flow of CF 4 gas. The (Helium) 
backside gas pressure can be adjusted and monitored using backside gas delivery system 26, 
which includes a pressure regulator. In addition, the substrate temperature can be monitored 
using temperature monitoring system 27. 

[0065] In an alternate embodiment, controllable process parameters can include a film 
material viscosity, a film material surface tension, an exposure time, a depth of focus, etc. 
[0066] With continuing reference to the table in FIG. 11, data can be recorded and stored 
digitally on controller 55 as a data matrix X , wherein each column in the matrix X 
corresponds to a given variation in a controllable process parameter (column in the table of 
FIG. 11) and each row in the matrix X corresponds to a specific spatial component. Hence, 
a matrix X assembled from the data in FIG. 11 has the dimensions 16 by 12, or more 
generally, m by n. Once the data is stored in the matrix, the data can be mean-centered 
and/or normalized, if desired. Centering the data stored in a matrix column involves 
computing the mean value of the column elements and subtracting it from each element. 


WO 03/058699 PCT/US02/38989 

16 

Moreover, the data residing in a column of the matrix can be normalized by the standard 
deviation of the data in the column. The following sections will now discuss the methods by 
which one determines the extent to which variations in controllable process parameters 
contribute to the spectral signature of spatial components. 

[0067] In order to determine the inter-relationships between variations in controllable process 
parameters and the spatial components, the matrix X is subject to multivariate analysis. Li 
one embodiment, principal components analysis (PCA) is employed to derive a correlation 
structure within matrix X by approximating matrix X with a matrix product {TP T ) of lower 
dimensions plus an error matrix E , viz. 

[0068] X=TP T +E, W 

[0069] where f is a (m by p) matrix of scores that summarizes the X -variables and P is a 

(n by p, where p<n) matrix of loadings showing the influence of the variables. 

[0070] In general, the loadings matrix P can be shown to comprise the eigenvectors of the 

covariance matrix of X , where the covariance matrix S can be shown to be 

[0071] 5 = X T X. (2) 

[0072] The covariance matrix S is a real, symmetric matrix and, therefore, it can be 

described as 

[0073] S=UAU T , & 

[0074] where the real, symmetric eigenvector matrix U comprises the normalized 
eigenvectors as columns and A is a diagonal matrix comprising the eigenvalues 
corresponding to each eigenvector along the diagonal. Using equations (1) and (3) (for a full 
matrix of p=n; i.e. no error matrix), one can show that 

[0075] P = U (4) 

[0076] and 
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[0077] f T f = A. & 

[0078] A consequence of the above eigenanalysis is that each eigenvalue represents the 
variance of the data in the direction of the corresponding eigenvector within n-dimensional 
space. Hence, the largest eigenvalue corresponds to the greatest variance in the data within 
the n-dimensional space whereas the smallest eigenvalue represents the smallest variance in 
the data. By definition, all eigenvectors are orthogonal and, therefore, the second largest 
eigenvalue corresponds to the second greatest variance in the data in the direction of the 
corresponding eigenvector which is, of course, normal to the direction of the first 
eigenvector. In general, for such analysis, the first three to four largest eigenvalues are 
chosen to approximate the data and, as a result of the approximation, an error E is introduced 
to the representation in equation (1). In summary, once the set of eigenvalues and their 
corresponding eigenvectors are determined, a set of the largest eigenvalues can be chosen and 
the error matrix E of equation (1) can be determined. 

[0079] An example of commercially available software which supports PCA modeling is 
SIMCA-P 8.0; for further details, see the User's Manual (User Guide to SIMCA-P 8.0: A 
new standard in multivariate data analysis . Umetrics AB, Version 8.0, September 1999). The 
contents of the manual are incorporated herein by reference. Using SMCA-P 8.0, for 
example, with the data of FIG. 11, one can determine the scores matrix T and the loadings 
matrix P , as well as additional information regarding the ability of each component to 
describe each variable in X and the total variation of each variable in X by a component. 
FIG. 12 presents the cumulative sum of squares R2X (cum.) of all of the variables in X 
explained by the extracted principal component(s) for the first three principal components 
and the cumulative sum of the total variation of each variable in X that can be predicted by 
the extracted principal component(s) for the first three principal components. 
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[0080] FIG. 13A presents the scores for each spatial component in t(l), t(2) space provided 
the exemplary data of FIG. 11 and FIG. 13B presents the loadings for each variable in p(l), 
p(2) space provided the exemplary data of FIG. 11. The data of FIG. 13A, in t(l)-t(2) space, 
displays the data variability through a measure of dispersion from the data center where, in 
particular, spatial components 1 and 2 are shown to reside outside the Hotelling T2 (5%) 
ellipse. This result indicates one should investigate the first and second principal components 
as shown in FIG. 13B and should further consider also components 3 and 4. From FIG. 13B, 
one can derive that the variations in controllable process parameters that would lead to a 
reduction in the magnitude of the spatial components could potentially be increasing the 
cooling gas pressure (i.e. helium backside pressure), decreasing the substrate holder 
temperature, decreasing the process pressure, decreasing the RF power and utilizing a 20 mm 
focus ring. 

[0081] Moreover, FIG. 14A presents the scores for each spatial component in t(l), t(3) space 
provided the exemplary data of FIG. 11 and FIG. 14B presents the loadings for each variable 
in p(l), p(3) space provided the exemplary data of FIG. 11. A similar conclusion can be 
drawn from analysis of FIGs. 14A and 14B and, therefore, the results of this analysis for 
generating a reduction of the spatial components is summarized in the table of FIG. 15. 
[0082] Utilizing the multivariate analysis summarized in FIG. 15 in conjunction with the data 
of FIG. 1 1 , one can reduce the data set of FIG. 1 1 to a more manageable set of data shown in 
the table of FIG. 16A. From the table of FIG. 16A and the (baseline) signature presented in 
FIGs. 6A.B, FIG. 16B presents a measured signature (baseline condition) and a correction 
(subtracted condition) signature according to the multivariate analysis, and FIG. 16C presents 
the difference signature once the correction (subtracted) signature is removed from the 
measured signature. After adjusting the controllable process parameters following the 
guidelines of the multivariate analysis to affect the difference signature of FIG. 16C, an 
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improved spatial uniformity of the scan of data for a process performance parameter is 
achieved and shown relative to the nominally measured scan of data in FIG. 17. In FIG. 17, 
the uniformity is improved by more than an order of magnitude (i.e. approximately 5% to 
0.5%). 

[0083] In an alternate embodiment, the implementation of multivariate analysis to determine 
a relationship between the controllable process parameters and the spatial components of 
process performance parameters can be achieved via design of experiment (DOE) 
methodologies. DOE methodologies are well known to those skilled in the art of experiment 
design. 

[0084] With reference now to FIG. 18A, a method of characterizing a material processing 
system according to an embodiment of the present invention is presented. The method 500 is 
described as a flow chart beginning with step 510 in which a controllable process parameter 
associated with a process performed in the material processing system is varied. The process 
performed in the material processing system can be, for example, the act of processing a 
substrate using a material processing system such as, for example, one of those described in 
FIGs. 1 through 5. In step 520, a scan of data, the data comprising a process performance 
parameter (PPP) as discussed above (i.e. etch rate, deposition rate, etc.), at, for example, two 
or more points above the substrate is measured and recorded. In step 530, the scan of data is 
transformed into spectral space. In step 540, a characterization of the material processing 
system is performed by identifying a process signature of the process performance parameter 
using one or more spatial components. Thereafter, in step 550, the process signature can be 
recorded in a data matrix as, for example, a column in a data matrix. 

[0085] In step 560, a decision is made as to whether an additional controllable process 
parameter should be varied. In order to further characterize the material processing system, 
steps 510 through 540 can be repeated, wherein an additional controllable process parameter 
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associated with a process performed in the material processing system is varied, an additional 
scan of data comprising a measurement of a process performance parameter is measured, an 
additional number of spatial components is determined from transforming the additional scan 
of data, and the material processing system is re-characterized by including an additional 
process signature comprising an additional number of spatial components. Furthermore, as 
before, the process signature can be stored in an additional column of the matrix in step 550. 
[0086] In step 570, the data assembled in the data matrix can be further processed utilizing 
multivariate analysis in order to determine inter-relationships between the variations in the 
controllable process parameters and the spatial components. Examples of multivariate 
analysis are principal components analysis (PCA) and design of experiment (DOE), which 
are described above. 

[0087] Referring now to FIG. 18B, a method for optimizing a process in the material 
processing system is described. In the method, a reference signature, deemed optimum for 
the given process performed in the material processing system, can be obtained. Utilizing the 
inter-relationships between variations in controllable process parameters and spatial 
components, the process is tuned by adjusting at least one controllable process parameter in 
step 610. In steps 620, 630 and 640, a scan of data corresponding to a process performance 
parameter is measured (step 620), the scan of data is transformed into spectral space to form a 
number of spatial components (step 630), and the resulting process signature is verified (step 
640). In the verifying step 640, the process signature is assessed to determine if the 
optimization of the process signature was successful. For example, if the optimal process is a 
uniform process, then the optimized process signature should comprise minimal amplitudes 
for each of its spatial components. If the verifying step 640 indicates a successful 
optimization, then the multivariate analysis is not altered in step 650 and a reference 
signature for the process in the material processing system is obtained in step 660. If the 
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verifying step 640 indicates an unsuccessful optimization, then the multivariate analysis can 
be altered and a series of steps described in FIG. 18A can be re-executed. 
[0088] Referring now to FIG. 18C, a method 700 for improving a process in a material 
processing system is described. In step 710, a relationship between the process signature and 
a controllable process parameter is determined. The relationship can be determined, for 
example, using data inspection, or any of the multivariate analyses described above (i.e. PCA, 
DOE, etc.). In step 720, a decision is made as to whether to improve the process. The 
improvement can, for example, involve improving the uniformity of the process performance 
parameter. In such a case, it would be advantageous to alter the process in order to minimize 
the amplitudes of at least one spatial component in the process signature, or minimize a 
difference signature formed from the subtraction of the process signature (FIG. 18A) from the 
reference signature (FIG. 18B). If no improvement is deemed necessary, all process data 
including the process recipe and process signature is recorded in step 730. If an improvement 
is deemed necessary, then the process is improved using a variation in at least one 
controllable process parameter in step 740. In step 750, a decision is made as to whether the 
method should be terminated. If not, the next process (i.e. next substrate, next batch, etc.) 
can proceed. 

[0089] With reference now to FIG. 18D, a method according to an embodiment of the present 
invention is presented. The method 800 is described as a flow chart beginning with step 810 
wherein the measured process signature (FIG. 18A) is compared with a reference signature 
(FIG. 18B) and a difference signature is computed as the difference between the measured 
process signature and the reference signature. The reference signature can be an ideal 
signature determined for the process or a signature obtained for a prior substrate. For 
example, the reference signature could constitute a spectrum of spatial components ideal for 
process uniformity. If the difference signature exceeds a threshold, then a process fault is 
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detected and the operator can be notified. Upon detection of a process fault, an alarm can be 
triggered in step 820 and a message can be sent to an operator (either locally or remotely) in 
step 830. 

[0090] In one embodiment, the threshold is exceeded when the difference signature for any 
one spatial component exceeds a fraction of the reference value for the spatial component. 
For example, when the aforesaid fractional deviation exceeds 20%, then the operator is 
alerted to perform the possibility of a process fault. The fault alarm can be indicative of the 
need for preventive maintenance, i.e. chamber cleaning or consumable replacement, or a 
defective process. In fact, a hierarchy of operator notifications can be employed, for 
example, a 10% fractional deviation can lead to an operator warning, a 15% fractional 
deviation can lead to an operator request for preventative maintenance scheduling and a 20% 
fractional deviation can lead to emergency notification for immediate attention. (Likewise, a 
programmed deviation may even be programmed to cause an automatic shut down of the 
malfunctioning equipment.) The controller 55 of FIGs. 1 through 5 can be connected via an 
integrated circuit (IC) manufacturing plant intranet to a centralized plant server for 
notification of a fab operator and/or an equipment supplier server located in the fab through 
which notification may be sent to the equipment supplier off-site via the internet. 
[0091] In an alternate embodiment, the threshold is exceeded when the sum of the spatial 
components comprising the difference signature exceeds a fraction of the sum of the spatial 
components comprising the reference signature. For example, when the aforesaid fractional 
deviation exceeds 20%, then the operator is alerted to perform the possibility of a process 
fault. 

[0092] After a decision as to whether to send a warning in step 820 is made, a decision can 
be made as to whether to improve the process in step 840. If an improvement is not 
necessary, all data including process recipe, spatial components, etc. can be recorded in step 
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850. If an improvement is deemed necessary, then a relationship between the difference 
signature and at least one controllable process parameter can be determined in step 860 and 
an improvement can be made through adjustment of at least one controllable process 
parameter in step 870. In step 880, a decision as to whether to end the process is made. If 
not, the next process or a new process can be executed while referring back to FIG. 18A. 
[0093] In the embodiments described herein, a one-dimensional scan of data has been 
utilized to determine a set of spatial components. Li an alternate embodiment, the scan of 
data can be multi-dimensional such as, for example, at least a two-dimensional scan of data. 
[0094] Although only certain exemplary embodiments of this invention have been described 
in detail above, those skilled in the art will readily appreciate that many modifications are 
possible in the exemplary embodiments without materially departing from the novel 
teachings and advantages of this invention. Accordingly, all such modifications are intended 
to be included within the scope of this invention. 
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CLAIMS: 

1 . A method of material processing, the method comprising the steps of 

characterizing a process, said characterizing comprising identifying a signature of said 
process wherein said signature comprises at least one spatial component; 
optimizing said process, said optimizing comprising identifying a reference signature; 
comparing said signature of said process with said reference signature for said 
process, wherein said comparing comprises determining a difference signature, and 
determining a process fault by comparing said difference signature with a threshold, 
wherein said process fault occurs when said threshold is exceeded. 

2. The method according to Claim 1, wherein said performing a process comprises 
processing a substrate. 

3. The method according to Claim 2, wherein said substrate is at least one of a wafer or a 

liquid crystal display. 

4. The method according to Claim 1. wherein said at least one process performance 
parameter is at least one of etch rate, deposition rate, etch selectivity, etch feature 
anisotropy, etch feature critical dimension, film property, plasma density, ion energy, 
concentration of chemical specie, temperature, pressure, mask film thickness, and mask 
pattern critical dimension. 

5. The method according to Claim 1, wherein said plurality of spatial components are 

Fourier harmonics. 

6. The method according to Claim 1, wherein said determining said relationship between 
said signature and said set of controllable process parameters comprises a multivariate 
analysis. 

7. The method according to Claim 6, wherein said multivariate analysis comprises principal 
components analysis. 
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8. The method according to Claim 6, wherein said multivariate analysis comprises design of 
experiment. 

9. The method according to Claim 1, wherein said at least one controllable process 
parameter comprises at least one of process pressure, RF power, gas flow rate, cooling 
gas pressure, focus ring, electrode spacing, temperature, film material viscosity, film 
material surface tension, exposure intensity, and depth of focus. 

10. The method according to Claim 1, wherein said improvement comprises an improvement 
of a spatial uniformity of said scan of data. 

1 1. The method according to Claim 1, wherein said improvement comprises a minimization 
of at least one spatial component. 

12. The method according to Claim 1, wherein said scan of data is a multi-dimensional scan 

of data. 

13. The method according to Claim 12, wherein said multi-dimensional scan of data is a two- 
dimensional scan of data. 

14. A system for material processing, the system comprising 
process chamber, 

device for measuring and adjusting at least one controllable process parameter, 
device for measuring at least one process performance parameter, and 

controller, said controller capable of characterizing a process, said characterizing 
comprising identifying a signature of said process wherein said signature comprises at 
least one spatial component; optimizing said process, said optimizing comprising 
identifying a reference signature; comparing said signature of said process with said 
reference signature for said process, wherein said comparing comprises determining a 
difference signature, and determining a process fault by comparing said difference 
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signature with a threshold, wherein said process fault occurs when said threshold is 
exceeded. 

15. The system according to Claim 14, wherein said process chamber is an etch chamber. 

16. The system according to Claim 14, wherein said process chamber is a deposition chamber 
comprising at least one of chemical vapor deposition and physical vapor deposition. 

17. The system according to Claim 14, wherein said process chamber is a photoresist coating 
chamber. 

18. The system according to Claim 14, wherein said process chamber is a dielectric coating 
chamber comprising at least one of a spin-on-glass system and a spin-on-dielectric 
system. 

19. The system according to Claim 14, wherein said process chamber is a photoresist 
patterning chamber. 

20. The system according to Claim 19, wherein said photoresist patterning chamber is an 
ultraviolet lithography system. 

21 . The system according to Claim 14, wherein said process chamber is a rapid thermal 

processing chamber. 

22. The system according to Claim 14, wherein said process chamber is a batch diffusion 
furnace. 
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